RUSSIAN ACADEMY OF SCIENCES 
PETERSBURG NUCLEAR PHYSICS INSTITUTE 



Preprint 2362 
June 2000 



S. I. Manayenkov 



New Method for Data Treating in Polarization 
Measurements 



Gatchina — 2000 



Abstract 



Precise formulas are derived for the expected values < £ >, < rj > and variances ££ 2 , 8rf 
of random variables £, rj describing the spin asymmetry of some reaction when a back- 
ground process contribution is negligible and appreciable, respectively. The variances of 
£ and 7] are proved to be finite. This property differs from that of the Caushy distri- 
bution which has an infinite variance. It is shown that < £ > is equal to the physical 
asymmetry which allows to find the asymmetry from experimental data without study- 
ing the detector efficiency This is the base of the proposed method of data treating. 
Asymptotic formulas for < rj > and 5rj 2 are also derived for a total number of events 
tending to infinity for a finite value of the background to signal ratio. 



1 Introduction 

The study of the polarization phenomena is of great importance in modern physics 
since it gives a valuable information about the spin structure of interaction under in- 
vestigation. Typical examples of the polarization measurements to be discussed in the 
present paper are the measurements of spin-spin asymmetries in inclusive deep inelastic 
scattering of leptons on nucleons (which give the spin-dependent structure functions of 
nucleons) or spin-spin asymmetries in semi-inclusive processes. A study of a decaying 
particle polarization in a final state of some reaction can be performed by measuring spin 
asymmetries in an angular distribution of detected particles. To get a high statistical 
accuracy of measured polarization observables we should integrate the result over all 
kinematic variables of which the studied physical quantity is independent. To perform 
such an integration we have to know the detector efficiency. Modern detectors are very 
complicated devices and a problem to describe an efficiency of the particle registration 
and the detector acceptance is very non-trivial. Therefore approaches in which one does 
not need the study of the detector efficiency have a great advantage. 

Let us consider some examples which explain our approach. To study the polarization 
of the A hyperon we may investigate the angular distribution of the pions and protons 
in the weak decay 

A°^p + n-. (1) 

It is well known that the proton angular distribution in the hyperon rest system looks 
like 

W(S,P) = —(l + aScos9) = — [1 + a(SP)/P] (2) 

47T 4:71 

where 9 is the angle between the proton three-momentum P (P = \P\) and the hyperon 
polarization vector S (S — \S\), a = 0.642 ± 0.013 U is the well known constant of the 
weak decay ([I]). The mean value of the detected number of the protons is 

An(K, P, S) = N A E(K, P, R)W(S, P)An P An K (3) 

where is the total number of the decaying hyperons, AQp is the solid angle of the 
registered protons and AQk denotes the solid angle of the decaying hyperon momenta, 
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E{K, P, R) is the detector efficiency. The latter depends on the hyperon momentum K, a 
position of the A decay vertex R and the proton momentum P. We would like to remark 
that we shall not discuss in the present work the problem of smearing which can, in 
principle, be important for the detector efficiency. We see from (|3|) that we have to know 
E(K, P, R), depending on many variables, to obtain the angular distribution W(S,P) 
from the measured quantity An/ (AQ P AQ K ). To reduce the statistical uncertainty of 
the obtained polarization of A we are to fit the angular distribution W(S,P) within the 
total detector acceptance for the proton registration. 

If it is possible to inverse the direction of the A polarization we may study the 
physical asymmetry 

An + + An_ 

where 

An± = An(K, P, ±S) . (5) 
It is easy to see from @, @ and (f|) that 

C = aS cos 9 . (6) 

We see from (|6|) that C is proportional to the hyperon polarization and does not depend 
on the detector efficiency if the solid angles AQk and AQp are small enough. In this 
approach we do not need to know the detector efficiency but we have encountered two 
problems. The quantities of An + and An_ in (Q) and (|j) are the expected values of the 
observed numbers p and n of decays (0) for the positive (p) and negative (n) hyperon 
polarization, respectively. If we consider the random variable £ 

f = ^ P) 
p + n 

then its expected value < £ > can be quite different from the physical asymmetry C 
given by (|J) or by the equivalent relation 

c = <P>-< n > (8) 
< p > + < n > 

for finite (and even small) values of < p > and < n >. We denote as < p > and 
< n > in (H) the expected values of the random variables p and n, respectively. We 
may suspect that < £ >— > C if An_ =< n >^ oo, An + =< p >^ oo. The former 
problem mentioned above is whether the difference < £ > — C is large or not at finite 
values of < p > and < n >. It will be shown later that < £ >= C if one may neglect 
the background contribution. 

In the approach under consideration we make use of only a small fraction of exper- 
imental events which belong to the solid angle AQp. The latter problem consists in 
utilization of all the observed events to reduce the statistical uncertainty of the obtained 
polarization of A . This problem can be non-trivial if the probability density of £ has 
the same property as the Cauchy density. Indeed, the Caushy random variable is equal 
to the ratio of two random variables having the Gaussian probability densities. It is well 



2 



known (see for example [0,0) that the Caushy distribution has an infinite second mo- 
ment and the probability to observe values of the Cauchy random variable which deviate 
significantly from the expected value is large. Moreover, the sample mean of N Cauchy 
random variables has the same probability density as every random variable. This means 
that one does not reduce the statistical uncertainty of the measured quantity considering 
the sample mean. But the random variable £ is just the ratio of two random variables 
p — n and p + n (see (|7|)). We shall show that the variance <5£ 2 of £ is finite if we use the 
conditional probability for £ and the Poisson densities for the random variables p and 
n. But it is just the case for the real experimental numbers of events which are positive 
integers. This result allows to use the total number of events to reduce the statistical 
uncertainty of the obtained polarization of A . It is well known how to do this. Indeed, 
let us divide the total kinematic region of Qp into N bins AQ\, AQ2, AQn and define 
the random variables £j, Q in the jth bin by the relations 

6 = ^, 0» 

Pj + n, 

S = I , (io) 

bj = a cos9j (11) 

where pj and rij denote the numbers of A decay events for the positive and negative 
hyperon polarization, respectively. In (]TTj) 9j is a value of the angle between the A 
polarization and the proton momentum in jth bin. Let us choose M bins (M < N) from 
N bins and denote them as 1th, 2nd, Mth. Since the expected values of are equal 
to C given by © (we are to replace 6 with 9j for the jth bin) the random variables Q for 
every j and ( have the expected values < Q > and < ( > equal to the A polarization 
S where 

M 

C = E&0 (12) 
i=i 

and the coefficients (3j are positive numbers obeying the relation 

M 

= (13) 

In a general case bj are given by expressions other than (|TT|) but they represent some 
known functions. They are used to make the expected values of all Q equal to each 
other. It is well known that (3j can be chosen to minimize the variance 5( 2 of the 
random variable (. The optimal choice looks like 

1 M 1 -1 

u Sj m =l °S m 



and gives the final result for S( 2 



M 1 

i=i 
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Some subtle points concerning properties of ( will be discussed in the next section. 

The method discussed above is applicable also for the case when the detector efficiency 
is unstable for the time of data taking. The longitudinal double spin asymmetry in 
inclusive or semi-inclusive deep inelastic scattering of leptons off nucleons is defined by 
the relation 

A _ dcr++ - da + _ 
da ++ + do + _ 

where da ++ (da + _) denotes the differential cross section of the studied process when 
both helicities of colliding particles are positive (have different signs). The typical time 
of data taking for the modern experiments is about few months. If the detector efficiency 
and the efficiency of the luminosity monitor changes considerably during this time we 
cannot use the relations 

d °++ = , d °++ = (17) 

where E denotes as before the detector efficiency, n ++ , n + _ are observed numbers of 
events and L ++ , L + _ are integrated luminosities when both helicities have the same sign 
and opposite signs, respectively. If the direction of the target (or beam) polarization 
changes for example every run we may ignore all the instabilities and define the random 
variables Q 



2j-l 2j 

n+ + — n+_ 

2j-l . 2j~ 



for two neighbour runs with numbers (2j — 1) and (2j). Formula ([18]) corresponds to 
relations (||), (0) with bj = 1. The expected value < Q > for every Q is equal to the 
physical asymmetry All and hence we may apply formulas (|T2D, (|I4"D, ( |i~5| ) to minimize 
the statistical uncertainty of the measured asymmetry. When writing ([T8|) we supposed 
that the integrated luminosities are equal to each other for the neighbour runs. 



2 Expected value and variance of measured asym- 
metry 

Let us consider the random variable £ defined by (^) and let the random variables p and 
n have the Poisson probability densities 

W p (p) = - x e~ x , W n (n) = V —e-y (19) 
pi n\ 

with the expected values for p and n being equal to x and y, respectively. First we 
consider the case when the background contribution is negligible. Since x (y) is propor- 
tional to the cross section for the case when the product of the initial particle helicities 
is positive (negative) the physical asymmetry looks like 

C = . (20) 

x + y v ' 
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It follows from (0) that £ becomes meaningless for p + n = 0, hence we should exclude 
the term with p + n = in computing its expected value < £ >. Then the formula for 
< £ > reads 



x p y n p - n _ x _ y 



W+ p ££ 1 p\ n\p + n 
x ^ x p ~ 1 y n 1 



e 



E 



-e 



-x-y 



n 



W +P t^>! (p-l)\n\p + 

JL y ^jT^ l _ e -*-v (21) 

W +p ^ lP \(n-l)\p + n • 1 } 

Putting p + n — 1 = m in the sums in and applying the Newton binomial theorem 
one gets 

<£>= 

e -x-y oo 1 r m yfn—n yn m x p yin-p 

~W^]^ m + ir]^(m-n)\ri~ y ]^^\(m-p)\ 

e ~x-y oo f x _^y\m _ y i _ e -x-y 

{x ~ v) ^\ J? = ^+y~> w + ■ (22) 

Since we have excluded in fl22|) the contribution with p = n = we are to divide the sum 
over p and n by the probability of events with p + n > which is 

W+ = 1 - Wp(0)W n (0) = 1 - e~ x - y . (23) 

Combining (^) with (^3|) we get our final result for the first moment of £ 

<i>= x -^y. (24) 

x + y 



Comparison of ( P4|) with (|20|) shows that the expected value of the random variable £ 
coincides with the physical asymmetry C not only for high statistics experiments when 
x ^> 1, y ^> 1 but for any values of x and y. The obtained result is rather surprising. It 
allows to measure the asymmetry in small kinematic regions (bins) and make use of total 
statistics for all the bins where the expected values of random variables are equal to each 
other, to improve statistical accuracy of the obtained asymmetry. As has been explained 
in the Introduction this method allows to avoid a study of the detector efficiency. To 
realize this program we need the formula for the variance of £ for every small bin to 
apply it in formula (|15|). 

To calculate the variance of £ let us consider the second moment < £ 2 > and represent 
it as the sum of four terms Ii, I2, I3 and I4 where 

< £2 >= / (p-n) 2 \ = / p(p-l) \ + / n(n-l) 
\ (p + n) 2 ' \(p + n) 2 ' \ (p + n) 2 

2 ((^) + (^)) = /l + /2 + /s + / - < 25 > 
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For calculation of I\ we make use of the relation 

1 1 



(p + n)i (q-l)\Jo 
for q = 2. This gives for I\ 



h 



-a(P+n) a Q-l da 



(26) 



(p + n) 



EE 



W+a&Jo (p-2)\n\ 



OO rr>P ni' 1 

U e - a{p+n ) adae -x-y 



x 2 e -x-y _oo m 



w. 



EE 



+ m=0 i 



X 



^ Jo (m-n)\n\ 



e~ a{m+2) ada 



-x-y 



poo 

/ ae~ 2a exp{(x + y)e~ a }da 
Jo 



(27) 



where W + has been defined in ( p3[ ) and m = p + n — 2. After integrating by parts integral 

[1 - e~ x ~ y - 



(|27|) can be represented as 

h 



x 



W + (x + y) 
where we denote by (f>(z) the function 

e* - 1 



(x + y)\ 



(2? 



/■Z p l _ ^ 

/ £ idt = e -*V 

Jo t ±^ 



/o t m • m! 

The integrals J2 and J3 can be calculated in an analogous way and we get 



(29) 



h + h + h 



(x - y)' 



{x + yf 



1 



(x + y) 



1 - e~ x -y 



(30) 



We have taken into account (|23| ) in (p0|). The calculation of ^4 is trivial. Indeed, remem- 
bering (|23|) and introducing m = p + n one has 

1 



p + n 



1 00 x p y n e~ x ~ v e~ x ~ y 00 1 m x m ~~ n y n 
W+ p J^> x p\ n\ p + n W+ ^im„? (m-n)!n! 



-x—y 00 



li 



E 



+ m=l 



(ar + y) T 



m • m! 



0Qz + y) 
1 - e- x -y 



(31) 



where 0(z) is defined by (p9|). Putting (^) and (31) into (|25|) and remembering (24) we 
get easily the final formula for the variance of £ 



{x + y) r (x-y) 2 



1 - e -x-y 
6 



(2 + y) 1 



(32) 



For numerical calculations with formula (|32|), it is convenient to use the integral repre- 
sentations for (f>(z) 

M z ) = - f Z ln(l - t/z)e- t dt , (33) 
Jo 

<f>(z) = l-e- z + z 2 [\l - t) ln(l - t)e~ zt dt . (34) 
Jo 

Now we are going to discuss some properties of the variance fl32|). For a high statistics 
case when i + ?/ > 1 we may use the asymptotic power series for (j)(z) valid if z — » oo 



£ tSt • (35) 

n=0 * 



Formula (|35|) can be easily derived at 2; 3> 1 from relation (|33|) through decomposition of 
ln(l— t/z) into a power series and integration over t with the limits and 00. Substituting 
( pop into (^) we get the asymptotic relation valid at x + y ^> 1 

^ 2 = ^- + 0((o: + y)- 2 ) (36) 
x + y 

where C is the physical asymmetry defined by fl20|) and 
0((x + y)~ 2 ) denotes terms of the order (x + y)~ 2 and smaller at x + y — ► 00. 
Since x and y are the expected values of observed numbers of events p and n, respec- 
tively, formula fl36|) says that <5£ decreases as reciprocal of a square root of the total 
event number. The solid line in Fig. 1 shows the dependence of <fi(z) on z, and the 
dash-dotted curve represents a dependence on 2 = x + y of the ratio of the variance 
calculated with fl32|) to its asymptotic expression given by d36|). We see that the ratio 



deviates significantly from the unity at z < 15 especially at z ~ 5. This means that 
widely used formula fl36"D for the variance of the measured asymmetry is not valid for 
the low statistics. We would like to stress that though £ is the ratio of the random 
variables p — n and p + n (see(|7|)) it has a finite statistical uncertainty at any x and y 
which follows immediately from (0), (|34]) and (|36"D . This result is in a contrast with the 
well known property of the Cauchy random variable having the infinite second moment. 
The reason for this is the used convention of removing the contribution of events with 
p + n = 0. Another difference of principle between the Cauchy distribution and the 
probability distribution of £ is as follows. Since for £ defined by (^) the numerator p — n 
and the denominator p + n are correlated we have for any integer numbers p, n and k a 
relation 

-l<(^) fc <l (37) 



•-p + TV 

which means that all moments of £ are finite if p+n > 0. The Cauchy random variable is a 
ratio of two uncorrelated random variables. We see also from (|32| ) and ( pDD that 5£ 2 = if 
|C| = 1. This property of 5£ 2 is due to the fact that for x = the probability to observe 
p > vanishes according to (|T9|), hence £ = — 1 for any n, which means <5£ 2 = 0. It 
follows also from (|l9|) that n = at 3/ = 0, hence £ = 1 for any p. The property under 
discussion can be proved from the obvious inequality 
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Fig. 1: Functions 4>(z), and x( z )- Solid and dashed curves show functions 4>(z), x( z ) 
defined by (|34|), (|36|), respectively. Dash-dotted curve represent the function z*(f)(z). 



1 >< £ 2 >=< [(£- < £ >)+ < £ >] 2 >= < f > 2 

which means that 

<5£ 2 < 1- < £ > 2 = 1 - C 2 . (38) 
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Let us consider the limit of expression fl32|) at x + y — > which will be used below. 
Applying ( j34|) at low x + y we have 

+ = x + y + 0{{x + yf) (39) 

and substituting (^) into (|3^) we get the relation valid at 
a; + ?/ — > 

^ 2 = l"(^) 2 = l-<e> 2 (40) 
K x + y' 

which realizes the upper limit in inequality (P5|). 

As has been discussed in the Introduction, to increase a statistical accuracy of the 
measured physical quantity (the polarization, the spin-spin asymmetry etc.) we are to 
utilize total number of experimental events. For this reason we consider all those bins 
in which the first moments of Q defined by @, fllPf ) are equal to each other. But in the 
practical application of formulas ([12]), (|14"D, (p"5|) we should remember that £j and Q are 
not defined for those bins where pj +n,j = 0. Hence we may use M bins among A" bins 
only (M < N). It is obvious that M cannot be larger than the observed number of the 
experimental events N exp . It is easy to get the upper limit for the variance of ( defined 
by (|12|) and (0). Let us consider very small bins so that A" ^> N exp . For this case the 
expected values of the observed numbers of events pj and rij are much less than unity 
(xj < 1, i/j < 1) and M coincides practically with N exp . We may apply fl4"U| ) for the 
variance of £j. Remembering the relation between £j and Q given by (|Hj) we get the 
following formula for the variance of Q valid if x 3 1, yj <C 1 

6$ = < > 2 ) = r^(l - h) < Q > 2 ) . (41) 

°3 3 

Taking into account that all < Q > are equal to < ( > we substitute (f4T|) into (|T5|) and 
come to the relation 

Since for physically interesting cases < < ( > 2 < 1 we get from ( f42|) the inequality 
of interest 

where < b 2 >m is the arithmetic mean of b 2 over M bins 

<b 2 >M=^J2b 2 . (44) 

Since M ~ N exp formula (f43|) shows that even for very small bins when observed event 
numbers pj, rij are about 1 we can get high statistical accuracy 
jV e ^p > 1. 

The procedure of extraction of the physical quantity < ( > and its variance 6( 2 from 
experimental data can be as follows. First, we calculate the arithmetic mean ( with the 
aid of the formula 

i M 

C = S E0 (45) 
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where Q has been defined by (^j) and (p~0|) . The arithmetic mean ( is an estimate of the 
physical quantity < ( >. After that we may estimate parameters Xj and yj putting 

Xj=Pj, Vj = nj (46) 
and apply fl4"6|) to calculate 5£? with the aid of relation ( j32|) 

^ ^a-^) ' 1 "^ 2 ' (47) 

which is an estimate of the variance of Q. Finally, we put 8(? into (JTSJ) instead of <5C? to 
get an estimate of the variance of ( where 

M _ 

^=[EV^]" • ( 48 ) 

Since p and n are integer numbers the random variable £ is defined for rational 
numbers in accordance with ([F]) . The probability to observe £ = (p — n) / (p + n) is 
not a smooth function of £. We may formally define the quantity f2 by the relation 

n(f , = ijy e (£ - - < < f + -) (49) 

which would become the probability density in the limit h — > if were a smooth 
function . Figure 2 shows the dependence of Q on £ at a; = 6, y = 2 for numbers of 
bins iV equal to 30 (top) and N = 150 (bottom) where h and N in ([49| ) are related with 
the formula h — 2/N. It is obvious that there is no smooth function Q given by fl4"9p 
at /i — > 0. Figure 3 shows that f2 corresponds to more smooth histograms than those 
presented in Fig. 2 if we increase x and ?/ keeping the asymmetry C (see (^)) equal 
to the same value as in Fig. 2. Figure 3 corresponds to higher statistics in comparison 
with Fig. 2 since x = 60, y = 20. It is easy to see from Fig. 3 that for high statistics Q 
resembles the Gaussian distribution. 

To find the confidence levels which correspond to the regions 

C - k5i < i < C + k8£ (50) 

{k = 1, 2, 3) usually used in the statistical analysis of the experimental data we are to 
calculate the cumulative distribution using the formula 

F(z) = W& <z)=T, X ^~ X ~ V /W + . (51) 



«<2 



pi n\ 



with W+ given by (|23|). The sum in fl51~|) runs over all positive integer p and n obeying 
the unequality z > £ = (p — ri) / (p + n) . The values of < £ >, 5£ 2 are calculated with 
the aid of @§, © and to compute w(JWf) (Jfe = 1, 2, ... ) 

= F(< £ > +^0 - F(< £ > -MO (52) 
relation ( ETj ) has been used. The results of the calculations are presented in Table 1. 
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0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 , 1 



Figure 2: Probability density Q for random variable £. Histograms are com- 
puted for x = 6 and y — 2. Bin number: top - iV = 30, bottom - N = 150. 



11 




0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 , 1 



Figure 3: Probability density Q for random variable £. Histograms are com- 
puted for x = 60 and y = 20. Bin number: top - N — 30, bottom - N = 150 
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The first line in Table 1 corresponds to the Gaussian distribution. We see that even 
at rather low statistics (x + y ~ 100) the numbers u(k5£) agree reasonably with the 
confidence levels for the normal density function. We see from the comparison of every 
line with the first line that the higher statistics the better is the agreement between the 
confidence levels for the random variable £ with those for the Gaussian distribution. It 
is easy to see from Table 1 that the agreement between the confidence levels becomes 
worse with an increase of | < £ > | (other things being equal) when the expected value 
goes close to its limits equal to ±1. 



X 


y 


<£> 


<5£ 




w(25£) 


w(3<5£) 










0.68268 


0. 


.94601 


0.99730 


3 


1 


0.5 


0.49712 


0.43339 


0. 


.95559 


0.96794 


12 


4 


0.5 


0.2242 


0.68386 


0. 


,95436 


0.99615 


120 


40 


0.5 


0.068682 


0.68344 


0. 


,95427 


0.99720 


7 


1 


0.75 


0.2535 


0.87089 


0. 


,96202 


0.99453 


28 


4 


0.75 


0.11886 


0.68360 


0. 


,95069 


0.99546 


140 


20 


0.75 


0.052456 


0.68165 


0. 


,95492 


0.99693 


19 


1 


0.9 


0.10016 


0.85166 


0. 


,95990 


0.99207 


76 


4 


0.9 


0.049045 


0.68116 


0. 


,94877 


0.99527 


380 


20 


0.9 


0.021822 


0.68228 


0. 


,95547 


0.99694 



3 Taking into account background 
events 

Usually measured numbers of events consist of events of the process under investigation 
and background events. Let the background events are described with random variables 
m, k both having the Poisson distributions 

W m (m) = —e- z , W k (k) = -eT* (53) 
ml k\ 

where m (k) describes the background contribution to the observed number p (n). Since 
the expected numbers of the events for the process under investigation are equal to x — z 
and y — t, then the formula for the physical asymmetry reads 

c = x - z - y + t . (54) 

x + y — z — t 

We suppose as before that the region of the kinematic variables is small enough therefore 
the detector efficiency can be considered as a constant and hence the ratio in ([54]) does 
not depend on the detector efficiency 
We can obtain the asymmetry if 

p + n — m — ky^O (55) 
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but we start our discussion considering the measured asymmetry in case when a total 
number of events p + n is larger than a total number of background events m + k 

p + n — m — k > . (56) 

The total probability to observe events in region ( p6| ) will be denoted as W + and is given 
by the relations 

W+=l-W-, (57) 



p n 771 j-k 

y „-x-y-z-t 



w V , , , , , 



q=0 



r=0q=0 



(z + t) q+r (x + y) 1 - 
(q + r)\ q\ 



e . 



(5f 



For transformation of the sums in (|58| ) we put q = p + n, I = m + k and after that we 
define r = I — q. We have applied in (|58|) the short notation 



e = e 

Making use of the well known formula |4[] 



-x—y—z—t 



E 

r=0 



r!(z/ + r) 



(59) 



(60) 



get from (]58 



for the Bessel functions I m (z) = (—i) m J m {iz) depending on the imaginary argument we 

(61) 



r=0 ' y 



-x—y—z—t 



Remembering another well known formula for the Bessel functions 



1 s , d« 



and putting 



u = CL\ 



ITpfexpW-!)} 



with arbitrary a obeying the inequality 



we get the useful expression 



a > 1 



(62) 



(63) 



(64) 



I n (2J(x + y)(z + t)) 



a n / Z + t \n/2 



(^_\ n/2 f n d<f>e- in +E(x + y,z + t, 
K x + y' Jo 



(65) 
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with 



E(x + y,z + t,<p) = exp{a(x + y)e i<t> + -(z + t)e~^} . 

n 



(66) 



Summing over r in (|6l|) with the aid of (|65|) and remembering ( 57|) we obtain the final 
result for W+ 



e f 27T 

1 / d(j)E(x + y, z + t, 

2n Jo 



1 - 



z + t 



-id> 



a(x + y) 

Let us define the random variable r\ by the relation 

(p — m) — (n — k) 
' (p — m) + [n — k) 

and consider its expected value in region fl56]) 

E/X p y n z m t k \ /p — n — m + k 
{—r-. — rrr (- r 
^ . , v p\ n\ ml kl ' v p + n — m — k 

p+n>m+k 1 1 



u p -u n -u m + u k , 



where for s = p, n, m, or k 

Us = e E 

p+n>m+k 

Making use of the obvious relation 

1 



x p y n z m t k 



(p + n — m — ky (j — 1)! Jo 
at j — 1 (valid in region (|56"D ) we get for U p 

f°° %P y n 2 m t 

f/„ = el da V 

P Jo _ , ±1 



n! m! &! ' + n — m — k 



p+n>m+k 



(p — 1)! n\ m\ k\ 



-a(p+n—m—k) 



xe \ dae 



,-yiS „ .71 yVTl -i-k 



s+n>m+fc 



s! n! m! fc! 



xe I dae~ a 

to 



a ^ (* + ?/) g £ (z +ty c _ a(Q _ l} = 

(x + y) l+r (z + t) 1 



q=0 I" Z=0 

oo oo 



r=Q «=0 



(Z + r)! Z! 



(67) 



(6* 



(69) 



(70) 



(71) 



(72) 



In ([7^) we put s = p — 1 and after that introduce q = s + n, I = m + k. We change the 
order of summation in (|72|) making use of r = g — I. We get the chain of the equalities 
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substituting (|60|) and then (|65|) into ( fT2| ) 

t^ = xe /°° rf«e- Q f] e" ar (^f ) § / r (2^/(x + + t)) 

• '() r=0 Z ~j~ X 

xe / dae~ a / V e- ra e~ ir ^ar r E{x + y,z + t } <j)) 
Jo Jo 2tt f^ 



xe 



2tt Jo 



2- 



da 



1 - e" a e- 



_E(x + y,z + t,< 



where E(x + y,z + t,<f>) has been defined by fl66|). The sum over r in 
due to (]5l|). Integrating over a in Q73D we get the relation 

c2tt 



u„ 



axe 



2tt Jo 



e^ln 



1 



1 - e-^/a 



E(x + y,z + t, 



(73) 

is convergent 
(74) 



The formula for U n can be obtained from ([74]) by the change x <-» y. The analogous 
calculation gives 



U, 



te r 2 ^ 



2na Jo 



"^ln 



1 - e-^/a 



E(x + y,z + t, 



(75) 



The relation for U m follows from ( [75] ) if we make use of the transformations t z. 
Putting ([FJ|), ( [75D and relations for Z7 n , t/ m into ( |69"D we get the final result for the 
expected value of the random variable rj valid if (|5T)|) is fulfilled 



< 77 >= 



2?r 



a(x - y)e i</> + 



t - 2 



ln ( 1 _ e -^ /a ) E ( :r + ^ Z + ^' 



(76) 



where W + is given by (|67]) , E(x + y, z + t,(j>) is defined by ( |66"D and e is given in (|59|). We 
would like to stress that < i] > does not depend on a if it obeys unequality (|64]). It is 
convenient for numerical calculations to put a close to the unity for a better convergency 
of integrals. To see that < rj > given by ( [fB]) does not coincide with the physical 
asymmetry C defined in for the case when the background process contribution is 
important, let us decompose ( [76]) into power series with respect to z and t. Considering 
them as small parameters and neglecting terms ~ t 2 , tz, z 2 and terms of higher orders 
it is easy to get the expression 



x — y . 

< 77 >= h 2(tx — yz) 

x + y 



(x + y) + e x y 



1 - e~ x ~y 



x + y 



+ ... . 



Expression Q77D can be compared with the power series for C 

C 



x — y (tx — yz) 

y + 2\ + ... 



(77) 



(7E 



x + y [x + y) 2 

which follows from ([54]). We see from a comparison of ( f77|) and ( |78|) that the difference 

r 0(x + y) + e - * - " 1 1 



< 77 > -C = 2(tx - yz) 



1 - e - x ~y 
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x + y (x + y) 2 



+ ... 



(79) 



is nonzero. Applying fl35|) we obtain from (Ifj3) the asymptotic expression for the differ- 
ence < 7] > —C valid at x + y 3> 1 



< i] > -C 



A(tx — yz) 
(x + y) 3 



+ 



10) 



Formula ( p0| ) becomes especially simple if there is no asymmetry for background events. 
Indeed, if z = t we get instead of (|80|) 

< n > c « 2 

C x + y x + ?/ 

Relation (^1|) shows that the fractional difference between the expected value of the 
random variable rj and the physical asymmetry C is proportional to two small factors: 
the background to signal ratio, 2t/(x + y) and the inverse value of the total number of 
events, (x + y) -1 . 

To calculate the variance of the random variable 77 defined by flBSj) we are to calculate 
< rj 2 > which (by definition) is given by the formula 



< 77 >=< 



p — m — n + 
p + n — m — k 



> 



W. 



£ ( 



x 



p y n z m t k \ (p — n — m + k\2 



p+n>m+k 



-)( 



p! n\ m\ k\ ' + n — m — k 



(82) 



We make use of formulas ([71]) with j = 2, (p]|), (^) to calculate sums over p, n, m, k 
in fl32|) as it has been demonstrated when obtaining the expected value of rj. The final 
formula for <rf > reads 



< Tj > 



2nW + Jo 



2k 



00 daae a E(x + y,z + t, 
1 - e- a e- ic t>/a 



(x - y) 2 e' a + (x + y) + —(z + t)e 



-2i<f> 



+ 



-(x-y)(t- z)e~^ + \(t- z) 2 e~ 3i(t> 



£ -^J^h(2^x + y){ z + t)) 
where W + , E(x + y, z + 1, <fi) and e are defined in fl67|) , (|66|) and ( |j9|) , respectively. To 



83) 



compute the variance of rj for region (|56|) we are to put (|83"D and (|76|) into the well known 
formula 

Srj 2 =< rj 2 > - < 77 > 2 . (84) 

Considering again z and t as small parameters, decomposing Srj 2 into power series 
with respect to z, t and retaining terms ~ z°, z 1 , t°, t 1 only we get for region fl56|) 



1 - 



( x ~ y ) 2 ] f ^ 
^x + v' ill 



+ 2/ 



+ (z + t) 



1 - 



- e~ x ~y 
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-(a; + 2 /)0(x + y)[l-(e^-l)- 2 ]- 



x + y 



} + 



+4 



(far - yz)(x - y) r 1 



x + y 



e x+y _ 1 

: y(x + -u) 

</>(x + + e~ x_y 



x + ?/ 



l - e ~ x -y 



15) 



where <f>(z) is given by ([34]) and the formula for reads 



o t 



m=l 



m z m! 



( \n 2 (l-t)e- zt dt. 
Jo 



(86) 



The behaviour of the function is presented in Fig. 1 with the dashed curve. To 
obtain the asymptotic behaviour of 5rf at large x + y, we are to use formula and 
the asymptotic power series for x( z ) at 2; ^ oo 



x(z) = £ 



m=2 



9m 



1 3 11 50 

z 2 z 3 z 4 z 5 



i7) 



where the coefficients g m in (|8y|) obey the recurrent relation 

g m = (m - l)g m -i + (m - 2)\ (88) 

with g 2 = 1. Substitution of fl3"5| ) and (^) into fl85|) leads to the relation of interest valid 
at x + y — > 00 



1 1 

+ 



g - y 

x + yJ ^x + y ' (x + y) 

, t + z ty 2 + zx 2 

+4 , . -4 . - +... . 



9) 



(x + y) 2 (x + y) A 

If we consider the expected value of t] in region (|55|) we are to make use of the relation 

<rj> (1-W Q ) = U P -U n -U m + U k + V p - V n -V n + V k , (90) 

where Wq denotes a total probability of events for which p + n = m + k. The formula 
for Wq follows immediately from ([H]) if we keep the term with r = in the sum and we 
have result 



W = £l (2J(x + y)(z + t) ) = ^- r d<pE{x + y,z + t 
K v ' Ztx Jo 



(91) 



In ( |9C|) U s denote quantities (s = p, n, m, k) which have been calculated above and V s 
can be defined by relation (|70| ) in which we are to sum over all p, n, m, k in the region 



p + n < m + k . 



(92) 
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It is easy to check that V p coincides with (— Uk) in which we are to make substitutions 
x «-> t, y <-* z. The expression for can be obtained from (—U p ) after the same 
substitutions. Applying the substitutions x <-> y, t <-> z one gets the formulas for V n 
and from the relations for V p and Vk, respectively. Putting expressions for Wo, U s , V s 
(s = p, n, m, k) into fl9~0|) we get the final formula for the expected value of r\ in region 
551) 



2tt Jo 



2- 



where e is defined in 



a(x - y)e i<s> + 
x 



ait 



In 



z)e 
1 



+ 



< 

a 

-y - 



-ic)> 



< T) > = 

E(x + y,z + t,4>) 
E(z + t,x + y,4>)} 



1 - e-^/c 
and + t, x + y, 



e/ (2V(a: + y)(2 + t)) 



is equal to 



(93) 



E(z + t,x + y, 



exp{a(t + z)e*+ + -(x + y)e~^} 
a 



in accordance with 

Decomposing (|93|) into power series with respect to z and t we get 



2(tx — j/z) - 



1 - e~ x -y 



[x + y) + e~ x ~ y 



<T]>- 
1 



X 



y 



x + y 



x + y 1 

Remembering (|78|) we get easily for the difference rj — C the relation 

-C = 



2(tx — yz) 



(x + y) 



-x-y _ 



(x + y)- 



< 7] > 
1 



e -x-y 



+ ... 



(94) 



(95) 



instead of (ff9|) . Formula ( Pq) shows that < i] > does not coincide with the physical 
asymmetry C in region (|55"D if the background contribution is not negligible. This is true 
in spite of applying the subtraction procedure: we consider p — m and n — k instead of 
p and n in fl68|). It is easy to check that the first term in the asymptotic formula for 
< rj > —C obtained from fl95"D coincides with fl80D . If one puts x, y, z, t equal to the 
observed numbers of the experimental events x = p, y = n, z = m, t = k, then formulas 
(|79|) and (|95|) can be used to estimate the systematic uncertainty due to the difference 
between the expected value of the random variable 77 and the physical asymmetry C. 
They can be improved if we make use of the precise expression for C and < 77 > given 

respectively. 



by (|54|) and 

To compute the second moment of rj in region ( |55| ) we may use formula ( [52]) if we add 
in the sum in the right hand side a contribution of events with p + n < m + k and put 
Wo instead of W-. This contribution is described by formula (|83| ) in which we should 
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use the substitution x *-*t, y 
like 



z, and W + — > 1 — Wq. The expression for < rf > looks 



2tt(1 - Wo) 



/. H (r 



-a p —i<j> 



< 7f > = 
£7(x + -2 + t, ( 



^ - y) 2 e- a + (x + y) + — (z + t)e 



+-{x -y){t- z)e- l<t> + — (t - zfe 



-3i(j) 



-2i</> 



+ 



2tt( 



£ r2n roo daae 01 \ 

1 ttft / ^ / (l _ a + 2, x + y, 4>) 

1 — Wq)Jo Jo v 1 — e a e ^/a' 



(t - zfe~ a + (t + z) + — (x + y)e 



+-(x - y)(t - z)e i4> + — (x - yfe 



-3i(f> 



-2i(j) 



+ 



w 



h{2J(x + y)(z + t) (x-y) 



z + t 
\] x + y 
Ix + y 



(96) 



with Wq defined in (pip. Putting (£56J) and (B^) into ( ]8~3| ) we come to the formula of 
interest for the variance of 77 in the region (|55|) . A decomposition of the formula for the 
variance into power series up to terms z and t gives the relation 



5 V 2 



/ x-y y-i f <, 



x + y) (^ + *) 



*x + y 

(x + y)(/)(x + y) 



+ 



— e - ^"^ ' 1 — e~ x ~ y 
1 - 2e- x ~y 



1 - e-^-f 



_ 4>(x + y) 
l - e- x ~y 

(x + y)e-*-y]} 



4(tx + yz) 
+ ~ Z7^7X{x + y) 



e -x-y 



4(tx — yz){x — y) \ 1 



x + y 



- <p(x +y)- e' x ~ v 



(97) 



(x + y)(l - e- x ~y) 

with <p(z) and x( z ) defined in fl34|) and fl86|), respectively. The asymptotics of (|9~?D at 
x + y ^> 1 is given by relation (|89|) . 

The most important limit of the obtained formulas corresponds to the case when 
x + y + z + tgoto infinity but the ratio Q — (z + t)/(x + y) is less than the unity. For 
this case the contribution of the region p + n < m + k is exponentially small. Indeed, 
making use of the well known asymptotic formula for the Bessel functions |4| 

' (98) 



\ / 2nz 

valid at z — > oo we get from (|6T| ) the relation 

Ty) hz + v)(z + t)]V*. 



w. 



m=0 



Z + tsm/2r e -x-y-z-t 



X 
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which shows that W- 



ex V {2^f(x + y)(z + t)} = 
yjx + y exp{-( v /x + y - sj z + t) 2 } 
47 (V* + y - >ftTz) [{x + y)(z + t)} 1 / 4 

exponentially if 

(y/xTy - Vz + t) 2 > 1 . 



(99) 



(100) 



Ignoring such exponentially small corrections we may obtain the asymptotic formula for 
< rj > which reads 



< rj > 



'X 



y + t 



x + y — t — z 
3(x + y + z + t) 2 



1 + 



(x + y - z - t) 4 
x — y — t + z 



+ 



[x 



y 



x + y + z + t 

+ ... 



(x + y 



(x + y — t — z) 



l 3(x + y + z + t) 
(x + y — z — t) 2 



:ior 



Comparison of ( p. 01] ) with (Q) shows that < t] > coincides with the physical asymmetry 

;t y+z+t tv « 1 • < io2 > 

[x + y — z — t) z 

It is easy to see that both conditions (|100|) and ( |102j ) are valid iix + y + z + t^>l and 
Q = (z + t)/(x + y) < 1. If conditions ( |100[ ) and ( |102| ) are fulfilled the variance of the 
random variable 77 looks like 



6 V 2 



x + y + z + t 
(x + y — z — ty 



1 + 



x 



X 



y-z + tf 



[x + y — z — ty 



yf 



ty 



[X 



(103) 



We have retained in ( |103| ) the greatest terms in the asymptotic power series only. For- 
mulas ( |101| ) and ( |103| ) correspond to the Gaussian distribution for the random variable 
77. In particular, relation ( |103| ) can be obtained from the widely used formula 



5 V 2 



2 = V + (^YSn 2 + ($LY8m 2 + (^) 2 6k 2 



\dpJ 



^dm' 



(104) 



where the dependence of 77 on p, n, m, k is given by fl68|) and dp 2 = p, 5n 2 = n, 5m 2 = m, 
5k 2 = k in accordance with the Poisson distributions fll9|), (|53|) for the random variables 
p, n, m, k. If one can neglect the background contribution (z — > 0, t — > 0) formula ( |103| ) 
reduces to (36). Formulas ( |101|) and (|103|) are obtained in Appendix. 

Figures 4 and 5 show the dependence of the ratio / = 
< rj > fC on A where C is defined by (154]) and the relation between x, y, z, t 
and A looks like 

x = x X , y = y X , z = z \ , t = t X (105) 
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the parameters x , y , z 0} t being chosen in such a way that x + yo + z + 1 = 1. The 
solid lines are calculated in region (|56| ) with the aid of formula ([F6|) and dashed curves 
are obtained for region ( |55|) using relation (p3[). The dash-dotted curves are computed 
with the asymptotic formula ( |101|) . Figure 4 shows that all the ratios go to the unity 
with an increase of the total statistics but the deviation from the unity increases if the 
mean value of the background to signal ratio b/s=(z+t)/(x+y-z-t) increases that follows 
from a comparison of Fig. 4a, Fig. 4b and Fig. 4c for which b/s = 12.5%, 33%, 75%, 
respectively. We see also that the higher the background to signal ratio the slower a value 
of < f] > jC tends to the unity with increasing A. The deviation from the unity for the 
curves presented in Fig. 4d is larger than for those of Fig. 4c though the background to 
signal ratio is the same for these two cases. It is due to the fact that the background 
asymmetry is equal to zero (z = t) for the curves shown in Fig. 4c while it has the 
opposite sign than the measured asymmetry (z — t < 0, C > 0) for the case presented 
in Fig. 4d. The behaviour of the ratio / =< rj > jC at small C (C ~ 1%) is presented 
in Fig. 5. Curves in Figs. 5a, 5b, 5c correspond to the physical asymmetry C = 0.025 
and b/s = 12.5% but the background events have zero, negative and positive asymmetry, 
respectively. We see that the greatest deviation of the curves from the unity corresponds 
to the case when C and the background asymmetry have opposite signs. The curves in 
Fig. 5d are calculated for the large background contribution {b/ s = 75%) nevertheless the 
maximum deviation of the curves in this case is even smaller than in Figs. 5b, 5c. This 
example shows that the asymmetry of the low background may lead to higher difference 
between < rj > and C than a high background without the asymmetry. This is true for 
low statistics only. We see from Figs. 4 and 5 that at large A the deviation of / from 
the unity is greater for the high background than for the low background. It is easy to 
see that the deviations of all curves from the unity is very appreciable. Figures 4 and 
5 illustrate a usefulness of formulas (ff6|), (|93|) which can be applied for estimating the 
difference between < rj > and the physical asymmetry C. We see from Figs. 4 and 5 
that the difference between formulas ([F(J) and fl93|) reaches about few per cents when 
X — x + y + z + t is about few dozens. It is easy to see also from a comparison of the 
curves presented in Figs. 4 and 5 that the asymptotic relation ( |101|) predicts < rj > with 
rather high accuracy for A > 20 — 30 even for the essential contribution of the background 
events when the deviation of < rj > from C is appreciable (~ 20%). 



4 Conclusions 

We have obtained the precise formulas for the expected values and the variances of 
the random variables £ and rj which correspond to the asymmetry of some reaction 
when the background process contribution is negligible and appreciable, respectively. 
Introducing the conditional probability distribution we have established that the 
expected value of £ is equal to the physical asymmetry C (see (0)) if there is no 
background process contribution to the reaction under study. We have obtained 
precise formula fl3"2|) for the variance of the random variables £ which is valid if one 



may ignore the background contribution. Relation fl3~2|) shows that the variance is fi- 



nite, and hence £ has nothing common with the property of the random variable with the 
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Fig. 4: Dependence of ratios < r\ > jC on A. Solid, dashed and dash-dotted curves 
are calculated using (|76|), ( p3| ) and (101), respectively. Parameters are: a) xq = 0.55, 
y = 0.35, z = 0.05, f = 0.05, C = 0.25, b/s = 0.125 b) x = 0.5, y = 0.3, z = 0.1, 
t = 0.1, C = 1/3, b/s = 1/3 c) x = 0.45, y = 0.25, z = 0.15, t = 0.15, C = 0.5, 
b/s = 0.75 d) x = 0.45, y = 0.25, z = 0.1, t = 0.2, C = 0.75, 6/s = 0.75. Parameters 
xq-, Uoi zq, to and A are defined in ( 105 ), C is given by (|54]) and 6/ 
denotes the background to signal ratio. 



(z + t)/(x + y-z-t) 
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Fig. 5: Dependence of ratios < r\ > fC on A. Solid, dashed and dash-dotted curves 
are calculated using (76), (93) and ( |101| ), respectively. Parameters are: a) xq = 0.46, 
y = 0.44, z = 0.05, t = 0.05, C = 0.025, b/s = 0.125 b) x = 0.45, y = 0.45, z = 0.04, 
t = 0.06, C = 0.025, b/s = 0.125 c) x = 0.47, y = 0.43, z = 0.06, t = 0.04, C = 0.025, 
b/s = 0.125 d) x = 0.355, y = 0.345, z = 0.15, t = 0.15, C = 0.025, 6/s = 0.75. All 
notations are the same as in Fig. 4. 



Cauchy density function. This is the result of defining £ in the region p + n > and 
making use of the conditional probability distribution. Formula (^) predicts that the 
variance vanishes if the physical asymmetry squared tends to the unity. For a large 
number of experimental events N exp the standard deviation <5£ goes to zero as \Jl/ N ex P 
according to (|5rj|). The finiteness of the variance of £ allows to utilize all the experimental 
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statistics accumulated in many bins (for which expected values of the random variables 
describing the measured asymmetry are equal to each other) to 

get a higher statistical accuracy of the asymmetry. We may consider the random 
variable ( defined by (|T2|) which has the expected value equal to the studied asymmetry 
but with the variance smaller than the variance 5Q which corresponds to the statistical 
uncertainty of the asymmetry obtained from the jth bin. 

As follows from ( ]75| ) and (^) the expected value of 77 both in region (^) and fl5"5| ) 
does not coincide with the physical asymmetry C (defined in (|54|)) for the process under 
investigation if the background contribution is not negligible. Relations ( [75] ) and Q§5D give 
the precise expressions for the expected values of the random variable 77 obtained with 
the conditional probability distributions defined in regions ( |56"|) and ([551), respectively. 



Exact formulas fl84|), fl83|) and (|76|) give the expression for the variance of 77 in region 
(|56|) . Precise relations (@), (|96"D and ( |93|) are to be used for computing the variance of 
the random variable 77 defined in region fl55|). We conclude from these relations that the 
variance is finite in both cases under discussion. This important result means that we 
may use the statistics of the experimental events accumulated in many bins to reduce 
the statistical uncertainty of the measured asymmetry. We would like to stress that, 
nevertheless, the obtained asymmetry has a systematic uncertainty since the expected 
value of rj deviates from the true physical asymmetry C in every bin. Formulas ([Ft]) and 
(|95f) can be applied to estimate the systematic errors due to the background contribution. 
The asymptotic formulas ( |101| ) and ( |103| ) for < 77 > and Srj 2 have been obtained for 
x + y + z + t ^> 1 and {z + t)/(x + y) < 1. This high statistics limit shows that 
< 77 > coincides with the physical asymmetry C and the standard deviation 5i] ~ 
y/x + y + z + t/(x + 2/-2-t)«lif conditions fllOOD and flID%) are fulfilled. 
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5 Appendix 

Let us consider the asymptotic behaviour of < 77 > and 5rf when x+y+z+t — ► 00 but the 
ratio Q = (z + t)/(x + y) < 1 . As has been mentioned the final results do not depend on 
the value of a but it is convenient to consider the limit of (|STD and ([76|) at a — > 1 (a > 1). 
Since the integrands in flB7p and ( [75] ) at <fi = 9 and <fi = 2ir — 9 are complex conjugated 
functions we may integrate from to it considering the real part of the integrals only. 
It is easy to see that the dominant contribution to the integrals in fl67|) and (|7q) comes 
from the regions < < (3 and 2ir — (3 < < 2tt with j3 2 ~ (x + y + z + 1)^ 1 1 if 
x + y + z + t is very large. Decomposing e ±lr ^ into power series up to terms <fi 2 we get 

1 i(t)Q 2 Q(1 + Q)- 



W+ » 1 - - r d<f)Re\ 

71 JO l 



Q (1-Q) 2 2(1 -Q) 



6 2 

exp{(x + y + z + t){\-^-) + i<f){x + y-z- t)}j. (106) 
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Neglecting exponentially small corrections it is possible to integrate in ( |106|) within the 
limits and oo. Making use of the following asymptotic formulas valid for positive and 
large u and v 



f°° 1 
K (u,v) — / exp{ — ut 2 + ivt}dt 
Jo 2 



u 3u 2 



V V v° 



f°° ] 

Ki(u,v) — / it exp{ — -ut 2 + ivt}dt 
Jo 2 



3u 



u 



V 



Ko(u,v) 



1 2% r 
(it) 2 exp{ — ut 2 + ivt}dt = — 

2 v 6 



l + 0( 



u 



(107) 



(which can be easily derived by integration by parts) and putting ( p.07|) into (|106|) we 
get that for large x + y + z + t 

W + = 1 (108) 

since all the integrals in ( 107 ) have no real parts. The corrections to formula ( 108|) are 
exponentially small in accordance with fl9"9"|). 

Since a behaviour of ln(l — e~ l< ^) at small looks like 



ln(l - e -*) = i- + ln0 - z| - ^ + ... (109) 

we need the asymptotics at large positive u and v for integrals containing logarithms 

f°° 1 
Gq(u,v) — / lntexp{ — ut 2 + ivt}dt = 
Jo 2 
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-m + - 3 m + ^(5) 
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K 1 (u,v) + 0(-) 
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/ (it) 2 hit exp{ — -wt 2 + ivt}dt 
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K 2 (u,v 
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mv 
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v° 



(110) 



with ■0(' z ) = r'(^)/r(z) being the digamma function. We shall obtain the expression 
for Gq(u,v). Other formulas ( |110| ) can be derived in an analogous way. To calculate 
Gq(u,v) we put the Frullani formula 
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into the definition of Gq in ( |110| ) and integrate over t 



dp 



Go(u,v) = / — / exp{ — -ut + ivt}e 



o (3 Jo 



— exp{ — -ut 2 + i(v + i/3)t} 



dt 
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One may write the integral 

M n = 

as the sum of the two integrals 
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iu 
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3iu 
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(112) 



dp 
o ~J 



(i + ip/vy 



dP f 
o ~P 
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r°° dp r 1 



i 



o P L(l + /3) n (l + ip/v) n 



(113) 



The first integral in (|113|) is equal to the digamma function due to the well known 
Dirichlet formula To calculate the second integral in ( |113| ) we write the chain of the 
equations 



r°° dp r 1 
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Jo 

(_l)n-l grv-l 


pt(l + p) n (l + i/3/v) n i 

f°°dPr 1 1 -I 


(n-l)! dx n - 1 J 


o /3 + /3 x + ip/yi x=1 
f°° dp \ i/v 1 i 


(n - 1)! ftr™- 1 J 
(-If- 1 d r 
(n — 1)! <9x 


o x i-x + ip/v x + pi x=1 

— - — lnf-1 =z lnt>. 

n-i Lx \ v Ji x =i 2 



(114) 



Putting the expressions for both integrals in ( |113] ) we get 



dP i 
o ~p 



1 



(l + z/3/^) r 
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?/>(n) + z— — lnt> 



(115) 



Combining ( |115|) with ( |112| ) we get the expression ( |110|) for G (u,v). 



To compute < rj >, we put (|108|) and (|109|) in (|76|) , decompose e ±l ^ into power series, 
integrating within limits and oo and making use of ( |107| ), ( |110| ) and get 
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{{x-y + t-z)(l- <p 2 /2) + z0(a; -y + z-t)} 

exp{ — — {x + y + z + t)<p 2 + i<p(x + y — z — ~ 

. .rlu 3u 2 1 i 

{x-y + t-z){- + — + — + — } 
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Formula ( |116| ) coincides with ( |101| ) since M = x + y + z + t and t> = x + y — z — t. We 
would like to note that the contributions to the real part give the terms i7i/2 and In (f> in 
the square brackets in ( |116[ ), their contributions being equal to each other in accordance 
with (jOTTD and fliTO!) . 

As has been mentioned above the contribution of region 
p + n < m + k vanishes in the limit x + y + z + t— > oo, Q — (z + t)/(x + y) < 1. This 
means that we are to get the same result (|116|) from (|93l) . Indeed, e/o in (p^) vanishes 
exponentially 



el {2yj{x + y){z + t)) 



exp{-x -y-z-t + 2J(x + y)(z + t)} -> . 



(117) 



A comparison of (|76|) and (^) shows that the additional terms in (|93|) after the substi- 
tution — > —0 (we remind that a = 1) look like (relation (|117 ) is taken into account) 



-Re 



[(x - y)e i * + (t- z)e- i(t> } 
E(x + y,z + t,(f>) ln(l - e irt> )d<p . 
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A comparison of ( |118| ) and ( |76|) shows that the only difference between them is the 
difference between the logarithms. The logarithm ln(l — e 1 ^) at small <fi looks like 



ln(l - e 



i(j>\ 



7T 



-i — h In 6 + % h 

2 ^ 2 24 



(119) 



As has been explained above only the terms m/2 and ln0 in the decomposition of the 
logarithm ln(l — e 1 ^) contribute to the real part of the integrals in (|118|) . But they 



have opposite signs in ( |119| ) (compare ( |119| ) with ( |109| )) and hence their contributions 
cancel each other totally. This means that (|93|) does really coincide with (^) in the limit 
x + y + z + t— > oo, Q < 1 (which means conditions ( |10(J| ) and ( |102[ ) to be valid). 

To find an asymptotic behaviour of < rf > at large x + y + z + t and Q = (z + 
t)/(x + y) < 1 we put relation ( |108[ ) in (^) and consider the limit a — » 1 (a > 1). Like 
in calculating < rj > the dominant contributions to the integrals in (|83"D come from the 



2s 



two regions < < j3, 2tt — (3 < < 2tt ((3 <C 1). These contributions are complex 
conjugate quantities. Hence we may consider the real part of the integral over in the 
limits and it. Since the integrals are convergent very rapidly we may use the limits 
and oo neglecting exponentially small corrections. To fulfill this program we have to 
establish the behaviour of the function 



Li 



1 - 



-da 



at small where 9 = e ^ in ( |120|) . The other functions of 
through L\. Indeed, it is obvious that 
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(120) 

can be expressed 
(121) 



Applying the new variable t 



e a and integrating several times by parts we get 
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It is easy to see that at small 
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-0 2 )(ln0 + ^) + ... 
1 4 V ' 4 
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-Hln0 + i-] + ... 



(123) 



and hence the derivatives dLf/d(j), d 2 L1/d<f) 2 , d 2 L\/d<f) 2 do not exist at = since 
expressions ( |123|) contain ln0. One may check that dL\/d(j), d 2 Ll/d<j) 2 are finite at 
= which means that LI may contain terms 0" In with n > 3 only. Differentiating 
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L\ over it can be easily obtained that 



6 3 Jo 



t 3 



-~ + C(2) + i<t>[-\ + C(2)] + $L[-\ + C(2)] + (124) 



where £(2) = 7r 2 /6 and denotes the Riemann (^-function. Putting ( |123| ) and 
into ( |122j ) we get 

L 1 = C(2) + #[C(2)-l] + ^[C(2) + i] 

+ (i(f)-(j) 2 )(\n(f) + i^) + ... . (125) 



The term eI\(2J (x + y)(z + t)) in fl33| ) is exponentially small owing to (0). Decompos- 
ing e ±%m<i> in (|83|) into power series, putting them and (|125|) , ( |121|) , (|108j) into (p3|) and 
integrating over <ft from to 00 one gets 

9 . s9 , 1 3m. m 3m 2 

< ^ >= x _ y + t _ z 2 + + + 

t> z t> 4 tr t> 4 
-^-^[(a;-,) 2 -(t-,) 2 ] (126) 

with -u = x + ?/ + ,2 + t and u = x + y — z — t. Substitution of ( [L26|) and ( |101| ) into ([54]) 
gives the final formula (|103|) for the variance 5rj 2 if conditions ( |100|) and ( |102|) are valid. 
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